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1. INTRODUCTION 

Time-delay is an inherent part of all practical control systems [1]-[3] and it may be caused by the time 
required for processing and/or transmitting of signals in the control loop [4], [5]. To design controllers for these 
group of systems, full understanding of the effects of time-delay on the system performance is required. Time- 
delay can reduce system quality of performance or cause system instability in worse case [6]. PID controller is 
generally used to control this class of systems due to its popularity [7], [8], simplicity [9], robustness and easy 
to use [10]-[12]. 

In PID controller design for time-delay systems, the first step is to establish the stability boundary for 
the system in its parameter space [13]. This is because such design is normally carried out using stability locus 
method [14]. This area of knowledge has been extensively researched. The method of computing all the 
stabilizing PID controller gains for a linear arbitrary order system with time-delay in the k;-kg plane with fixed 
kp was reported in [7]. More so, the procedure for the computation of the entire stability gains in the k,-k, 
plane with fixed kg was presented in [15]. The stabilization controller parameters for a time-delay integral 
fractional-order system under the control of fractional-order PID controller were determined in the k;-kp 
stability region for a given kg, as reported in [16]. The controller gains in the k,-k, plane for a given kg yielded 
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general stability region while the global stability region was obtained after sweeping over the permissible valus 
of ka. In [17], [18] parametric methods were used to establish the stability region in the k;-kg and kp-k; planes 
for a fixed k, and k4, respectively. Then GA was used to determine the optimum controller gains in the 
established region of the plane. 

PID parameters for second-order plant with time-delay system was obtained by [15] in the k,,-k; plane 
with a fixed kg for specified gain and phase margins using dominant poles method. PID controller tunning 
method for integrating system with time-delay and inverse response was developed in kp-k; stability region 
for a fixed kg in [19]. After determining the stability region, weighted geometrical center (WGC) approach 
was used to determine the required k, and k; gains for a given kg. Though the method produced good results, 
it is computationally intensive. PI-PD controller for time-delay systems was tuned by [20] using WGC method 
which yields satisfactory performance compared to some other methods in the literature. In [14] CCSR method 
for PI-PD controller design was proposed for unstable systems with time-delay. Both the experimental and 
simulation results shows that the method is superior to some of the others reported in the literature. The method 
is simple, less computational and with satisfactory performance. 

All the reviewed PID controller design methods were carried out either in the k;-kg plane with a fixed 
kp [21] or in the k,-k; plane for a fixed kg [22]. The CCSR method has several merits [14]; however, it is 
seldom applied in PID controller design. Its application in PI-PD controller design can be extended to PID 
controller design in the kp-k; plane for a fixed kg due to their k;-k, stability locus similarity. The method is 
amenable to the stability region in the k;-kq plane for fixed kp being a convex polygon. Despite the amenability 
of CCSR design method, its application has not been extensively explored in the research community, hence 
the need for this study. 

A PID controller design for time-delay system in k;-kq plane with a fixed kp using the CCSR method 
is reported in this study. The equations relating controller gains, system parameters and time-delay were 
derived; and then used for plotting the stability boundaries in k;-kg plane for fixed kp. Thereafter the centroids 
of the stability region, the required controller gains, were calculated. Using the calculated controller gain, the 
fixed k,,, the system and time-delay transfer functions, the closed-loop unit step response were plotted for all 
the generated stability regions on the same axes. The required TDPMs for characterising the designed systems 
were obtained from its closed loop step response plots. 

Based on the TDPMs obtained for the various combinations of controller gains, a unique gain combination 
can be selected for a system considering its performance specifications. In practice, conflicts do exist among the 
TDPMs; therefore, a compromise is made when selecting any of the controller gains combinations. 


2. METHOD 
2.1. Determination of stability region in the k;-kq plane for fixed kp 

The closed-loop transfer function required for the analysis and design of PID controlled time-delay 
system can be derived from the block diagram of unity feedback control system shown in Figure 1. In 
Figure 1, R(s), E(s) and Y(s) is the reference input, error, and output respectively. The expressions for the 
plant, time-delay and controller transfer function are given in (1), (2) and (3) respectively. 


N(s) 


G,(s) = 38 (1) 
where, N(s) and D(s) is the plant transfer function numerator and denominator respectively. 
Ga(s) =e™ (2) 
where, T is the time-delay in seconds. 
2 
G(s) = ad Sas ead) (3) 


Ss 


The expression in (1) was decomposed into its even and odd parts after substituting jw for s as shown 
in (4) for easy application of D-decomposition method [23]. The (—w?) term has been removed from 
N,(—w?), N,(—w?), D.(—w”) and D, (~w?) terms in (4) for compactness purpose. 


Net+jwWNo 
Det+jwDo 


GpGw) = (4) 
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Y(s) 


Gq(s) 


Figure 1. Block diagram of unity feedback time-delay control system 


The formulation of the stability boundaries in kp, k; and kg space is a three-dimensional problem. To 
simplify the controller design and analysis, it can be reduced to two-dimensions by fixing one of the gains and 
then finding the stability region in the plane of the remaining two gains. The overall stability region can then be 
determined using the stability regions in the plane of the two parameters by sweeping over the fixed parameter 
values. In the proposed method, k, was fixed and the stability boundary in k;-kg plane was determined. The 
equations and conditions required for determining the stability region in the k;-ką plane for a fixed k, are 
presented as follows [17], [18], [23], [24]: 


For w = 0 
k, =0 (5) 
For w > 0 
= (w?NoDotNeDe) cos(wT)+w(NoDe—NeDo) sin(wT) (6) 
p -(Ne?+w2No?) 
k, = w? (NoDe—NeDo) cos(wT)-w(NeDe+w?NoDo) sin(wT)+ki(Ne?+w?No?) (7) 
a w2(Ne2+w2No7) 
The frequencies w = wm, (M = 1,2,...... ) are the frequencies at which the line of a given value of k, 


intercepts the graph of k, versus w using (6). The stability boundary in the k;-kg plane is formed by the line 
obtained from (5) and the lines generated using (7) when w,, is substituted for w. It should be noted that m is 
the number of points of intersection of k, line with the graph of k, versus w, that is, the number of lines 
obtainable from (7). 


2.2. Determination of controller gains using CCSR method 

The advantages of this method over others are as follows: i) no iterative optimization process is required, 
ii) system’s closed loop stability is guaranteed, iii) trial and error are eliminated, iv) no error manipulation, v) low 
computational load, vi) controller can be design using only the stabilizing controller parameters region, vii) it can 
be used for any system with closed stability region and viii) it has high precision. 

Detailed discussion on the formulation and advantages of CCSR method of controller design can be 
found in [14]. This method is amenable to PID controller design for time-delay system in the k;-kq plane for 
fixed k, because, the stability region in this plane is a convex polygon (triangle in this study). 

To use this method, the coordinates of corner points of the stability boundary were determined. 
Assuming the number of corner points on the stability boundary is n, then the coordinates of the corner points 
in k;-kg plane can be represented by (kj1, kar), (kiz, Kaz), «0... (kin, Kan). The coordinate of the CCSR (kpc, 
kic), for fixed kp, which is the required controller gains in k;-kg plane can be calculated using (8) and (9) 
respectively [14]. 


n 
— W=tkpj 


kpc a (8) 


Xj=1Kij 
ki = (9) 
2.3. System performance 

There are two major states embedded in the time domain performance of feedback control systems, 


namely: i) transient-state, which describes the speed of the system closed-loop response to the step input and 
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ii) steady-state, which gives the accuracy of the system response under decayed transient condition. It should 
be noted that the system transient-state depends on swiftness and closeness by which the system responds to 
and tracks respectively the step input. The swiftness depends on the value of rise time (T,) and time-to-peak 
(Tp), while closeness depends on the value of settling time (Ts) and percentage overshoot (%0S) [25]. The 


lower the value of these TDPMs the better the system performance. The quality of system steady-state response 
depends on the value of steady-state error (ess) and the lower it is, the better. 


3. DEMONSTRATING EXAMPLES 


To demonstrate the proposed methods, three Examples were used. The detailed steps of PID controller 
design, testing and characterization using the proposed method are fully explained and presented for Example 
1 in section 3.1. To avoid repetition, these steps are skipped for Examples 2 and 3 as presented in sections 3.2 


and 3.3 respectively. 
3.1. Example 1 

Design of PID controller for an integrating second-order time-delay system with system transfer 
function given by (10) [26] was considered in this example. 


1 = 
G(s) = ae 3 (10) 


For this system, using (2), (4) and (10), Ne = 1, N, = 0, De = —w?, and D, = 1. 

The range of kg and k, for stability was first determined by plotting the trajectory of kg against kp 
for w=[0, œ] and k; = 0 as shown in Figure 2 using (6) and (7). Also based on (6), the plots of k, against œw 
are shown in Figures 3 for the determination of the relevant frequencies for any fixed value of kp. 
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Figure 2. Plotting k_d against k_p for Example 1 Figure 3. Plot of k_p against œ for Example 1 


It can be seen from Figures 2 and 3 that the range of k, is 0-1.717. To find the controller gains in the 
ki-ką plane the fixed gains considered based on range of k, as shown in Figure 3 are k,=0.2, 0.4, 0.6, 0.8, 1.0, 
1.2, 1.4 and 1.6. The detail design for k,,=0.2 is presented as follows: 

As earlier explained in subsection 2.1, when w=0, the stability boundary equation that can be used for 
generating one of the boundary lines was derived using (5). To generate the squations that can be used to form 
the remaining boundary lines, the frequency at the point of intersection of the line of k,=0.2 (with the graph 
of k, versus w in Figure 3) was determined. From Figure 3, the points of intersection are two and their 
corresponding frequencies are 0.322 and 1.991 rad/s, respectively. Substituting these frequency values into (7) 


yielded the two equations required for generating the remaining stability boundary lines. The resulting stability 
boundary equations are shown in (11)—(13). 


k, =0 (11) 
kaı = 9.6447k; — 0.8467 (12) 
kaz = 0.2523k; + 2.2257 (13) 
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Based on these equations and the similar equations obtained when the remaining fixed k, values were 
used following the steps used for k,=0.2, the system stability boundaries in the k;-kg planes shown in 
Figures 4 and 5 was obtained for k,=0.2 and for the complete set of fixed kp considered respectively. 


Convex stability region for kp = 0.2 


kd 
kd 


—©— Stability region boundary 


cA WE (kit, kd1)=(0, 2.226) 05 
-1F He (ki2, kd2)=(0.331, 2.346) F 
Á (ki3, kd3)=(0, -0.8467) 
1.5 1 1 l 1 f 1 fl i 4 
-0.05 o 0.05 0.1 0.15 0.2 0.25 0.3 0.35 3 
ki ki 


Figure 4. The stability boundary in the k_i-k_d Figure 5. The stability boundary in the k;-kg plane 
plane for k_p=0.2 for Example 1 for all set of fixed k, for Example 1 


In this example, the explanation given in subsection 2.2 was used for determining the controller gains 
as the CCSR in the k;-kg plane. From Figure 4, the stability boundary is a triangle, therefore there are 3 corners 
making n=3. The coordinates of the corner points as indicated in Figure 4 are (kj1, kg )=(0, 2.226), (kiz, 
ka2)=(0.331, 2.346) and (kiz, kqg3)=(0, -0.8467). Based on these coordinate points, using (8) and (9) the 
coordinate of the CCSR which is the required gains was determined to be (kic, Kg-)=(0.1103, 1.2418). The 
location of these controller gain values is indicated hexagram marker in Figure 6. 

The obtained kic and kac and the corresponding fixed kp were used in conjunction with the plant and 
time-delay transfer functions to plot the system closed loop unit step response. The required TDPMs for the 
characterization of the designed system were obtained from the step response graph. The adopted TDPMs used 
in this study are rise time (T,), peak time (Tp), percentage overshoot (%OS), percentage undershoot (%US), 
settling time (T) and steady-state error (ess). 


3.2. Example 2 
An integrating second-order time-delay system with left hand side zero with transfer function shown 
in (14) [19] was treated in this example. 


0.6(—0.3s+1)e— 925 
s(s+1) 


G(s) = (14) 
For this system Ne = 0.6, Nọ = —0.18, De = —w?, and D, = 1. 
The PID controller design in the k;-kg plane was considered for the following fixed k, values:1, 2, 3, 
4, 5, and 6. The stability boundary for these range of fixed k, is shown in Figure 7. The obtained kic and kac 
in the stability region associated with each of the fixed k, were also used accordingly as explained for Example 
1 in section 3.1. The closed loop unit step response graph was plotted and the necessary TDPMs obtained. 


3.3. Example 3 
In this example a second-order time-delay system with left hand side zero whose transfer function is 
shown in (15) [23] was considered in this example. 


G(s) = (-0.5st1)e~ 9-65 (15) 


(s+1)(2s+1) 


For this system N, = 1, N, = —0.5, De = 1 — 2w?, and D, = 3. 
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In this example, the PID controller design in the k;-kg plane was considered for the following fixed 
kp values: kp=-0.5, 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, and 3.5. By applying the method used for Example 1 the 
stability boundary for the considered range of fixed k, shown in Figure 8 was obtained. The obtained kiç and 
kac in the stability region associated with each of the fixed k, were appropriately used as earlier explained. 
The closed loop unit step response graph was plotted and the needed TDPMs obtained. 
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Figure 6. The CCSR in the k;-kg plane for k,=0.2 Figure 7. The stability boundary in the k;-kq plane 
for Example 1 for the set of fixed k, for Example 2 
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Figure 8. The stability boundary in the k;-kg plane for the set of fixed k, for Example 3 


4. RESULTS AND DISCUSSION 

The fixed k, values and the corresponding kic and kg, gains in the k;j-kg plane for Examples 1, 2, 
and 3 are presented in Tables 1(a), 1(b), and 1(c) respectively, where it can be seen that k;, and kg, increase 
with the value of k, in all the three (3) examples considered. These changes in the controller gains are a form 
of controller tunning, because for every change a new combination of controller gains is generated in the 
stability region. For better understanding of the effects of changes in controller gains on the system response, 
closed-loop unit step response graphs for the controlled and uncontrolled system (UCS) were plotted using the 
same axes for each of the three illustrated examples. 

The system closed-loop unit step responses in Examples 1, 2, and 3 are shown in Figures 9, 10, and 11 
respectively. As shown in Figures 9 and 10 (Examples 1 and 2), the response of the UCS is better than the response 
of some of the controlled systems. Therefore, the required TDPMs value for both the controlled and UCSs were 
determined for performance analysis. On the other hand, in Figure 11, the response of the UCS can never attain the 
reference input, hence its TDPMs value cannot be used for performance analysis. The TDPMs for each controller 
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gains combination using the fixed gains and UCS (where applicable) as indicator are presented in the bar charts 
shown in Figures 12(a) and 12(b), 13(a) and 13(b), and 14(a) and 14(b) for Examples 1, 2, and 3 respectively. 


Table 1. Controller gains in k;-kg plane convex stability regions for (a) Example 1, (b) Example 2, and 
(c) Example 3 


(a) 


kp 0.2 0.4 0.6 0.8 1.0 1.2 1.4 


kic 0.110 0.219 0.330 0.440 0.555 0.667 0.775 
kac 1.229 1.284 1.345 1.396 1.457 1.509 1.558 


(b) 


kp ! 2 3 4 5 6 


ki, 0.833 1.667 2.520 3.333 4.167 5 
kac 2.640 2.799 2.972 3.117 3.276 3.432 


© 


kp -0.5 0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 


kic 0.208 0.416 0.625 0.833 1.046 1.258 1.473 1.689 1.918 
kac 1.070 1.233 1.403 1.571 1.746 1.916 2.087 2.256 2.432 


Response 


4 1 A 1 1 1 4 1 1 
o 2 4 6 8 10 12 14 16 18 20 
Time (seconds) 


Figure 9. Unit step response graph for Example 1 


Response 


4 a L 1 L 4. L L 1 1 
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Figure 10. Unit step response graph for Example 2 


It can be observed from Figure 12(a), that the minimum and maximum values of %OS corresponds to 
fixed k, of 0.2 and 1.6 respectively and the minimum and maximum values of T, correspond to fixed k, of 1.2 
and 0.2 respectively. Also, from Figure 12(b) the minimum value of T, corresponds to fixed kp of 1.6 while 
the maximum corresponds to fixed k, of 0.2 and 0.4 but the minimum and maximum value of T, correspond 
to fixed kp of 0.6 and 0.2 respectively. From Figure 12(a) the only scenario with %US was when k,=1.6. In 
Figure 12(b) the system of Example 1 yielded good steady-state response for the designed controllers as ess=0 
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except for when fixed k,=0.2 and for the UCS that gives ess of 0.016 and 0.004 respectively which can be 
tolerated in practice. 


Response 


o 5 10 15 20 
Time (seconds) 


Figure 11. Unit step response graph for Example 3 
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14 i 


12 


TDPMs Value 
n 


Figure 12. Bar chart of TDPMs for fixed k, for Example 1, in (a) TDPMs are %0S, %US and T; and in 


(b) TDPMs are T,., T, and ess 


From Figure 13(a) the minimum and maximum values of %0S were attained when the system was 
uncontrolled and at fixed k, of 6 respectively while the minimum and maximum values of T; correspond to 
fixed k, of 5 and 1 respectively. From Figure 13(b), the minimum value of T, corresponds to fixed kp of 6 and 
the maximum was attained for UCS. But the minimum and maximum values of T, correspond to fixed k, of 3 
and 1 respectively. Also revealed in Figure 13(a) is the system %US which has its minimum and maximum 
values corresponding to UCS and fixed k, of 6 respectively. It can be observed from Figure 13(b) that the 
steady-state response for the system of Example 2 is equally good because ess=0 for the designed systems 
except for when kp=1 that produced ess of 0.004 which will have negligible practical implications. 

From Figure 14(a), the minimum and maximum values of %OS correspond to fixed kp of 1.5 and 3.5 
respectively and the minimum and maximum values of T, correspond to fixed kp of 2.5 and -0.5 respectively. 
Also, from Figure 14(b) the minimum value of T, corresponds to fixed kp of 3.5 while the maximum value 
was at fixed k, of -0.5 and 0. But the minimum and maximum values of T, corresponds to fixed k, of 3.5 and 
-0.5 respectively. In Figure 14(a), the minimum and maximum values of %US was obtained when the fixed 
kpwas -0.5 and 3.5 respectively. From Figure 14(b), the steady-state response of the system of Example 3 is 
better because the designed systems give ess of 0 except for the system with fixed k, of -0.5 that gives ess of 
0.01 which can still be accommodated from practical point of view. 

A situation where a system has the lowest value for both the %08S and T, was not encountered in the 
results. Also, simultaneous low values of T, and Tp were not realized in any of designed systems. From these 
observations and explanation on the system performance in subsection 2.3, it can be inferred that, none of the 
systems can be said to have absolute good transient response since non is having absolute swift response to 
step input and absolute close tracking of the step input by the response. 
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Figure 13. Bar chart of TDPMs for fixed k, for Example 2, in (a) TDPMs are %0S, %US and T, and in (b) 
TDPMs are T,, T, and ess 
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Figure 14. Bar chart of TDPMs for fixed k, for Example 3, in (a) TDPMs are %0S, %US and T; and in (b) 
TDPMs are T,, T, and ess 


In control engineering practice, different applications normally require different system performance 
specifications in the form of TDPMs. Therefore, the selection of PID controller gains combination will depend 
on the required system performance specifications. The selection can only be made among these designed 
controllers because the controllers fall under the region where the systems stability is guaranteed. Most of the 
times, conflicts do exist among the TDPMs, therefore, to select any of the controller gains combinations, 
compromise must be made among any of the conflicting or competing TDPMs. It should be noted that, the 
number of controller gains combinations which determines the number of possible controller selection can be 
increased by increasing the number of fixed k, considered in the system stability region. 


5. CONCLUSION 

PID controller design procedure for time-delay systems using the CCSR method in the k;-kg plane 
for fixed value of kwas presented and more so, the system stability region in the said plane was established. 
In all the cases considered, the values of the determined controller gains (k;, and kac) increase with that of 
fixed kp. The performance analysis was based on system transient and steady-state which make it simple to 
understand. Based on the obtained TDPMs and system performance specifications, a unique combination of 
controller gains can be selected for the system to improve its performance. Again, selection of any of the 
controller gains combinations, normally requires a compromise among any of the conflicting TDPMs. The 
number of gains combinations determines the number of possible controller selections, and the former can be 
increased by increasing the number of fixed k, in the system stability region. 
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